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Abstract 

Using the field-theoretical methods we studied the evolution from BCS de- 
scription of a non- Fermi superconductor to that of Bose-Einstein condensation 
(BEC) in one loop approximation. We showed that the repulsive interaction 
between composite bosons is determined by the exponent a of the Ander- 
son propagator in a two dimensional model. For a / the crossover is also 
continous and for a = we obtain the case of the Fermi liquid. 
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INTRODUCTION 



The problem of the crossover from BCS superconducting state to a Bose-Einstein con- 
densate (BEC) of local pairsEHi becomed very important in the context of high temperature 
superconductors (HTSC). While at the present time there is no quantitative microscopic 
theory for the occurrence of the superconducting state in the doped antiferromagnetic ma- 
terials, it is generally accepted that the superconducting state can be described in therms 
of a pairing picture. The short coherence length (£ ~ 10-20 A) increased the interest for the 
problemtHil because it showed that the BCS equations of highly overlapping pairs, or the 
description in terms of composite bosons cannot describe the whole regime between weak and 
strong coupling. The mean field method developed by different authorsiHB and solved ana- 
litically in two and three dimension, and the Ginzburg-Landau descriptionBSiil showed that 
the evolution between the two limits is continous, no singularities during this evolution ap- 
pearing. The zero temperature coherence length in the framework of field-theoretical method 



has been in Ref. [12|. The problem of the BCS-BEC crossover in arbitrary dimension d using 
the field-theoretical method has been extesively discussed in Refs. |I7)|-|T8|, where the chemical 
potential, the number of condensed pairs and the repulsive interaction between pairs have 
been calculating using the analogy with the field-theoretical description of superfluidity. 

In this paper we apply this method to study the crossover problem for a non-Fermi 
superconductor described by the Anderson modelffl (See also Refs. IjHU), to study the 
crossover between weak coupling and strong coupling. The paper is organized as follows. 
In Section II we study the weak coupling model for a d — 2 non-Fermi superconductor. 
Section III contains the strong coupling limit. To make the paper self-contained we present 
in Appendix the Lagrangian formalism for the superfluid phase following Refs. [T~5|JT6| . The 
results are discussed in Section IV. 



WEAK COUPLING LIMIT 

The BCS-like model for the non-Fermi system is described by the Lagrangian 

c = tfGo v T + ^(g- 1 )*^ - ao^ t v;^t (i) 

where the normal state is described by the Green function 

where u> c < uj < u> c , g(a) = ira/ (2 sin (7ra/2)) and A < is the coupling constant describing 
the attraction between electrons. The Green function given by Eq. (0) contain a cut-off uj c 



2 



and the exponent a. This form has been proposed first by AndersonEl to describe the 2D 
non-Fermi properties of the superconducting state. 
If we introduce the two-component fermionic field 



(3) 



the non-interacting part of the Lagrangian ([!]) is 



: = *t(^ o \ 



(4) 



In order to calculate the partition function 

Z = J D^D^exp iJ^C 
we will transform the interaction contribution from the Lagrangian (|1]) as 



cxp 



(5) 
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where f x — J dt J d d x as in Ref. |T6|-|T8|, and A = X tj;^ is a bosonic field. The partition 
function defined by Eq. (||) will be expressed using Eq. (§) in a bilinear form as 



Z = D&DV / DA^DA exp — / A f A 



A, 



o J x 



exp 



-A 
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Performing the integral over the Grassmann fields the partition function becomes 

Z = J DA^DA exp (iS ef/ [A\A] + — f A 1 " a) 
where S e ff[A\ A] is the one loop effective action, which can be written as 

\ -A T /(a)(Po + UPj) 



(9) 



where f(a) = u"g (a) and the trace Tr has been used according to the meaning from Ref. 



I6|. 



In the mean field approximation the integral from Eq. (|8]) can be performed using the 
solution given by the saddle point and for T ^ the critical temperature T c will be obtained 

as 

1/2q 

(10) 
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where A(a) = g 2 {a)2 2a sin (tt(1 - «))/tt, C(a) = T 2 (a)[l - 2 1 ~ 2a ]C(a), .0(a) = T(l - 
2a)r(a)/(2ar(l — a)), T(x) being the Euler's gamma function. This expression is valid 
only in the limit < a < 0.5 and a positive critical temperature implies for the coupling 

(uj c /uD) a /{A(a)D(a)). We have to mention that 



constant the condition I A, 



> A c , with \ c 



the critical temperature obtained in Eq. (TTD) , calculated also in Ref. ^ is different from the 
one obtained in Ref. |22|p^ , |25| , and it is easy to show that it gives the exact BCS result in 
the limit a —>■ 0. If we consider the effective action as 

/(^(po-ccp)) 1 - o 

o /(oOCpo + ^p)) 1 -" 



S eff [A\ A] = -iTr In 



iTr In 



IAI 



/(^(pS-eW-*. 

and the system as space time independent the partition function can be written as 



Zq exp 



An 



-A f A 



(12) 



Z$ containing the non-interacting contribution, and we get for the renormalized coupling 
constant A the expression 



1 

A 



1 



P 



d 2 p r dpo 
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Using the integral 



we calculated A as 



where 



k0 (k 2 -E 2 + t V y 
i 

A = 



Ai 



47rag~ 
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Ao Ai 
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2vr (pt-eip)) 1 -* 

T(/-l/2) 1 



(13) 



r(Z) e 21 ~ 1 



2<\ 



(14) 



(15) 



cos (vr(« - 1)) 5(1/2, 1/2 - a) \uj d J 

B(x,y) being the Euler beta function B(x,y) = T(x)T(y)/T(x + y). The expression given 
by Eq. dl5|) is positive for a < 1/2. The new coupling constant A has to be also negative 
in order to have superconductivity (A < 0) and this condition is satisfied if |Ao| < A i . If 
we consider also the condition A c < |Ao| we get the general condition for the bare coupling 
constant A , A c < |A | < Ai, which is satisfied for < a < 0.5. 

We mention that for the weak coupling limit Aq — > 0~ the BCS limit studied in Ref. 



24] is reobtained, but we also showed that the critical constant calculated from the critical 
temperature is smaller than Ai(a), which also satisfies condition Xi(a — > 0) = 0. 

In the limit Ao — > — oo, called the strong coupling limit, we expect an important effect of 
the non- Fermi character of the electrons in the coupling constant. 
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STRONG COUPLING LIMIT 



In this limit we consider A(x) = A + A(x) and consider the action S e ff[AJ, A] obtained 
from Eq. (0) as 



where 



S eff [A\A] = -iTr In 



1 + GnA 



Go 1 



/(a)(po-^(p)) 
-At 



l-a 



/(«)(p-o+c(p)) 



1-Q 



(16) 



(17) 



A 



-A 
-At 



which can be written as 



oo -i 

S eff [A\A] = -iTrJ2~ GqA 
i=i 1 



with 



Go(po,p) 



-A 
-At 



P(a)(p 2 -e(p)y- a -W 
We are interested in quadratic terms in A and we will take the approximation 

S e// [At,A] = 5g ) / (0) + Sg ) / (q) 
which contains the quadratic contributions. The first term in Eq. (|2T|) has the form 



Si%(0) = ^iTr 



1 



2 f 2 (a)(p 2 - e(p)) 1 -" - |A| 2 



(A 2 AtAt + A t2 AA + 2|A| 2 |A| 



H — iTr- 



:2|A| 



2 /^(pg-e^p)) 1 -"-^! 2 ' 

which will be approximated, taking in the dominator A « as 



1 1 

o( 2 ) fn\ ~ _»7V 

efA) 2 ! r / 4 (a)(^-e 2 (p)) 2 ^ 



A 2 At At + At AA + 2|A| 2 |A| 2 



(19) 



(20) 

(21) 

(22) 
(23) 



(24) 

the last term giving no contribution to the renormalized coupling constant. Following the 
same approximation we calculated S^fJq) as 



-iTr A At 

2 f 2 (a)(p - aP)Y- a (Po + Qo + £(p + q)) 1 -" 



H — iTr 



2 / 2 (a)(p + C(p)) 1 - a (Po + go - e(p + q)) 1 "" 



(25) 
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From Eqs. ([24]) and fl25|) we have 

£W(0H _ B (V2.3/2 a -2 a ) (2m) 3^ 



X 



/(^ ^ + l)^ [ A2AtAt+At " M+2IAI2IAI; 



(2vr) 2 (p 2 + m£ a ) 3 " 



and 



£ (2) _ sin ( 7r ( 1 ~ a))B(a, a) ml _ 2a f d 2 p 1 



47r/ 2 (a) 7 (27r) 2 (p 2 + me a + q m + q 2 /A) 1 - 

sm(n(l — a))B(a,a) 1 _ 2a f d 2 p 1 



-2a 



47r/ 2 (a) 7 (27r) 2 (p 2 + me a - q m + q 2 jA) x - 



2a 



The integrals from Eqs. ([^) and (27) can be performed using the formula 



1 T(N-d/2) 1 



(p 2 + A 2 ) N (ATi) d / 2 T(N) (A 2 ) N ~ d / 2 



and we obtain 

C^= m 22 " 4 " B ^^-^) ["A 2 A t A t + A f2 AA + 2|A| 2 |A|' 

167r 2 / 2 (a) 1 - 2a £ 2 " 4a L ^ 

m sin (7r(a — 1)) B(a,a) AA* 

+ 167r 2 / 2 (a) 2^ (e a + g + g 2 /4m)- 2a 

m sin (7r(a — 1)) B(a,a) A^A 

+ 167r 2 / 2 (a) 2^ (e a -q Q + q 2 /Am)- 2a 

Using the approximation 



and using the notation 
we obtain from Eq. (pq) 



A 

At 



where 



£^ = 167r2 ^ 2(a) sin (tt(1 - a))B(a, a^^M* 



-/^o -go - ^ - Mo 
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nib = 2m being the boson mass and /j,q the chemical potential 

1 2 2 - 4Q E(l/2,3/2-2 Q ) 
^° / 2 (a) (1 -2a) sin (vr(l -«))£(«,«)' 1 a K ' 



The velocity c of the sound mode is 

.2-4aRfi / 9 3/9 _ 9^ 

Al 2 ^ - 1 (33) 



/10 1 2 2 - 4 ^(l/2,3/2-2a) ,,,,.,„_, 



° m b f 2 { a ) (1 — 2a) sin (tt(1 — a))B(a, a)m 
and the repulsive interaction Aob between pairs is 



7T 



2 



2 4 ~ 4a J B(l/2,3/2 - 2a) 



Ao6(a) m(l-2a)[sin(7r(l-a))5(a,a)] 2 



We mention that lim Q ^ ^ob(ct) = 2n/m a result identical to the result obtained in Ref. [L6 
for the two dimensional case. 



I. RESULTS AND DISCUSSIONS 



Using the field-theoretical methods we studied the crossover between BCS and BEC in a 
non-Fermi liquid. The weak coupling case leads to the same results as in the mean field like 
mo delsiHH In the strong coupling limit we showed that the pairs form a Bose gas with a 
repullsive coupling constant which is controled by a. 
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APPENDIX: 



In this appendix we briefly present the field-theoretical formulation of the Bogoliubov 
theory for the interacting Bose gas. The system of an interacting Bose gas is described by 
the Lagrangean 

£ = i {$* [ido - i + /i ] $ - A |$| 4 + c.c} (Al) 

where $ is the operator corresponding to the complex scalar field, e the kinetic energy 
operator, /i the chemical potential and A the bare repulsive coupling constant. The model 
described by Lagrangean (|A1|) has a global symmetry to the transformation 

$(x) -> e ia $(x) (A2) 
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where a is a constant defining the transformation. At T = this symmetry is spontaneously 
broken and it is associated with the occurrence of the superfluid phase. According to the 
Goldstone theorem the dynamical restoring of the symmetry implies the occurrence of the 
Golstone bosons and it can be included in the theory taking 

= e i6(x) ($ + $(£)) (A3) 

However, in the simple mode which is important for the pair description we will neglect the 
Goldstone mode and will take 6(x) = 0. The expression for the Lagrangean (|A1|) becomes 

£ = I [id - e + no] <£ + [-id - e + /i ] $*0) 

- A ($ 2 <T(x) 2 + 5> 2 (x)$ 2 ) - A (& &(x) + <S>*(x) 2 $ 2 )} (A4) 

In the p-representation using 




Eq. (|A4j) becomes 



where 



and 



£ = -$*(x)M$(x) (A6) 

u | Po-£(p)-^)-4Ao|$o| 2 -2A |$ | 2 , 

' ! -2A |$*| 2 - Po -s(p)-U(x)-AX \%\ 2 1 

C/(x) = d 6(x) + ^- [W(x)] 2 (A8) 

Eq. (|A7p has been obtained neglecting the terms containing V 2 #(x), which is irrelevant in the 
low momentum approximation, and the term jV#(x), where j is the current associated with 
x- variation of the phase. For the elementary excitations spectrum we neglect U (x) (which is 
equivalent to neglect the x-dependence of 6(x)) and from the condition 

detM = (A9) 

we get 

E 2 (p) = e 2 (p) + 2ti e(p) = e 2 (p) + 4A |$ | 2 ^(p) (A10) 
The spectrum expressed by ( |A10| ) can be approximated, in the limit p — > 0, as 

E = u \p\ (All) 



with uq = \J jjo/m, is gapless and it remains gapless to all orders in perturbation theory. For 
large momentum from ( |A10|) we get 

E(p) =e(p) + 2A |$o| 2 (A12) 
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